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Q ■ Abstract 



Photon diffractive dissociation, ^p — ► Xp, has been studied at HERA with the ZEUS 
Q^ ' detector using ep interactions where the virtuality Q 2 of the exchanged photon is smaller 

than 0.02 GeV 2 . The squared four-momentum t exchanged at the proton vertex was 
determined in the range 0.073 < \t\ < 0.40 GeV 2 by measuring the scattered proton in the 
ZEUS Leading Proton Spectrometer. In the photon-proton centre-of-mass energy interval 
176 < W < 225 GeV and for masses of the dissociated photon system 4 < Mx < 32 GeV, 
the t distribution has an exponential shape, dN/d\t\ oc exp (— b\t\), with a slope parameter 
b = 6.8 ± 0.9 (stat.) ±}% (syst.) GeV" 2 . 
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1 Introduction 

The reaction 7p — > Xp, in which the photon diffractively dissociates into an hadronic state 
X with mass Mx, has been investigated with real photons at a photon-proton centre-of-mass 
energy W of about 14 GeV Q. Recently it has also been studied at HERA using the process 
ep — > eXp for photon virtualities Q 2 < 0.02 GeV 2 and W ~ 200 GeV @, ||. The comparison 
of the fixed target data and the HERA data indicates that the dissociation of real photons 
has similar characteristics to the dissociation of hadrons, as expected in the framework of 
Vector Meson Dominance (VMD) j|, [5|. In this model, the photon is assumed to fluctuate 
into a virtual vector meson prior to the interaction with the proton. The interaction can be 
described by Regge phenomenology || and, at high energy, is dominated by the exchange 
of an object with the quantum numbers of the vacuum, referred to as the pomeron. An 
exponential fall of the differential cross section da/d\t\ ex exp (— b\t\), at small values of \t\, is a 
typical feature of diffraction; here t is the square of the four-momentum transfer at the proton 
vertex. Regge theory also predicts that the diffractive peak shrinks as W increases according 
to b = b + 2a' In (W 2 /M x ), where b and a' are constants || 0. 

The studies of diffractive real photon dissociation at HERA have so far focussed on the shape 
of the Mx spectrum |2|, |3|. The t distribution for the reaction jp — > Xp has been measured 
only by the fixed target experiment M, which found that in the range 1.4 < Mx < 3 GeV and 
for 0.02 < \t\ < 0.1 GeV 2 the t dependence is exponential with a t-slope 5^5 GeV -2 . At 
HERA, measurements of the t distribution have been performed for the diffractive dissociation 
of virtual photons in the range 5 < Q 2 < 20 GeV 2 ||, and for elastic vector meson production, 
7£> — > Vp, both with real and with virtual photons [9-18]. In all cases the distribution has 
an approximately exponential shape. The t-slope is b = 7.2 ± 1.1 (stat.) l a 'g (syst.) GeV~ 2 
for the diffractive dissociation of virtual photons at (Q 2 ) = 8 GeV 2 . For elastic p° production 
b depends only weakly on W but varies from approximately 10 GeV -2 for Q 2 w [9-11] to 
approximately 5-7 GeV~ 2 for (Q 2 ) = 10 GeV 2 [12], |13| . It is therefore interesting to extend 
these measurements to diffractive real photon dissociation. 

In this paper we report the first determination at HERA of the t distribution for the process 
7P — ^ Xp, where 7 is a photon with Q 2 < 0.02 GeV 2 . The present measurement is based on a 
sample of photoproduction events collected using the reaction ep — > eXp at W ~ 200 GeV |3[]. 
The sample was defined by the requirement that the scattered positron be measured in a 
calorimeter close to the outgoing positron beam line and a final state proton carrying at least 
97% of the incoming proton momentum be detected in the ZEUS Leading Proton Spectrometer 
(LPS) [jTTj| . The LPS was also used to measure the transverse momentum of the proton, from 
which t was calculated. This is a technique similar to that used to measure the t distribution in 
the photoproduction of p° mesons, 7^ — > p°p [11], and for the diffractive dissociation of virtual 
photons, 7*p — ► Xp M. 



2 Experimental set-up 
2.1 HERA 

The data presented here were collected in 1994 at HERA which operated with 820 GeV protons 
and 27.5 GeV positrons. The proton and positron beams each contained 153 colliding bunches, 

1 



together with 17 additional unpaired proton and 15 unpaired positron bunches. These addi- 
tional bunches were used for background studies. The integrated luminosity for the present 
study, which required the LPS to be in operation, is 0.9 pb _1 . 



2.2 The ZEUS detector 

A detailed description of the ZEUS detector can be found elsewhere JT9|, EDI. A brief outline of 



the components which are most relevant for this analysis is given below. Throughout this paper 
the standard ZEUS coordinate system is used, which has the origin at the nominal interaction 
point, the Z axis pointing in the proton beam direction, hereafter referred to as "forward", the 
X axis pointing horizontally towards the centre of HERA and the Y axis pointing upwards. 
The polar angle 9 is defined with respect to the Z direction. 

Charged particles are measured by the inner tracking detectors which operate in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. Immediately surrounding the 
beam-pipe is the vertex detector (VXD), a drift chamber which consists of 120 radial cells, 
each with 12 sense wires [pi]] . It is surrounded in turn by the central tracking detector (CTD), 
which consists of 72 cylindrical drift chamber layers, organized into 9 superlayers covering the 



polar angle region 15° < 9 < 164° [22 



For the energy measurement the high resolution depleted-uranium scintillator calorimeter 
(CAL) is used [23]. It is divided into three parts, forward (FCAL) covering the pseudorapiditytl 



region 4.3 > r] > 1.1, barrel (BCAL) covering the central region 1.1 > rj > —0.75 and rear 
(RCAL) covering the backward region —0.75 > 77 > —3.8. Holes of 20 x 20 cm 2 in the centre of 
FCAL and RCAL accommodate the HERA beam-pipe. Each of the calorimeter parts is subdi- 
vided into towers, which in turn are segmented longitudinally into electromagnetic (EMC) and 
hadronic (HAC) sections. These sections are further subdivided into cells, which are read out by 
two photomultiplier tubes. Under test beam conditions, the energy resolution of the calorime- 



ter was measured to be oeJE = 0.18/JE(GeV) for electrons and (Je/E = 0.35/JE(GeV) for 
hadrons. The calorimeter noise, dominated by the uranium radioactivity, is in the range 15- 
19 MeV for an EMC cell and 24-30 MeV for a HAC cell. 

The luminosity is determined from the rate of the Bethe-Heitler process, ep — > e'yp, where 
the photon is measured with a lead-scintillator sandwich calorimeter (LUMI-7) located at 
Z = —107 m in the HERA tunnel downstream of the interaction point in the direction of 



the outgoing positrons ||23 |. A similar calorimeter (LUMI-e) at Z = —35 m detects positrons 
scattered at very small angles. In this analysis, the LUMI-e was used to tag photoproduction 
events with positrons scattered at angles up to about 5 mrad and to measure the scattered 
positron energy, E' e . The LUMI-e covers the range 7 < E' < 21 GeV. The energy resolution of 



both calorimeters is (JeJE = 0.18/JE(GeV). 



The Leading Proton Spectrometer (LPS) [JTTJ] detects charged particles scattered at small 
angles and carrying a substantial fraction, Xl, of the incoming proton momentum; these parti- 
cles remain in the beam-pipe and their trajectory is measured by a system of silicon micro-strip 
detectors very close (typically a few mm) to the proton beam. The detectors are grouped in six 
stations, SI to S6, placed along the beam line in the direction of the outgoing protons, at 23.8 m, 



x The pseudorapidity rj is defined as r/ — — In (tan (9/2)). 



40.3 in, 44.5 m, 63.0 m, 81.2 m and 90.0 m from the interaction point. The track deflections 
induced by the magnets in the proton beam line allow a momentum analysis of the scattered 
proton. For the present measurements, only the stations S4, S5 and S6 were used. With this 
configuration, for xl close to unity, resolutions of 0.4% on the longitudinal momentum and 
5 MeV on the transverse momentum have been achieved. The effective transverse momentum 
resolution is, however, dominated by the intrinsic transverse momentum spread of the proton 
beam at the interaction point which is ~ 40 MeV in the horizontal plane and ~ 90 MeV in 
the vertical plane. For xl close to unity, the LPS covers the range 0.25 < pt < 0.65 GeV, 
where pt is the transverse momentum of the proton with respect to the incoming beam direc- 
tion. As discussed previously fll ], the incoming beam direction and the beam position with 



respect to each station are determined using the reaction ep — > ep°p at Q 2 ~ 0. Protons with 
Pt < 0.2 GeV and xl ~ 1 are too close to the beam to be measured. For the events considered 
here the geometric acceptance of the LPS is approximately 6%. 



3 Data selection and background subtraction 

3.1 Trigger 

ZEUS uses a three-level trigger system |TJ| ^(J . At the first-level a coincidence between signals 
in the LUMI-e and in the RCAL was required. An energy deposit greater than 5 GeV was 
required in the LUMI-e. In the RCAL the deposit had to be larger than 464 MeV (excluding 
the towers immediately adjacent to the beam-pipe) or 1250 MeV (including those towers). 
The angular acceptance of the LUMI-e limits the Q 2 range to the region Q 2 < 0.02 GeV 2 . 
The small RCAL threshold essentially selects all photoproduction events. The second and 
third trigger levels were mainly used to reject beam related background. Parts of the data 
stream were prescaled [^, [25| in order to reduce the high event rate resulting from the large 
photoproduction cross section. 



3.2 Reconstruction of variables 

The photon-proton centre-of-mass energy squared, W 2 = (q + p) 2 , where q and p are the 
virtual photon and the proton four-momenta, respectively, was determined by W 2 ~ ys, with 
y « E^jE e = [E e — E' e )/E e and s the squared centre-of-mass energy of the positron-proton 
system; here E 1 is the energy of the exchanged photon and E e denotes the energy of the 
incoming positron. The W resolution is 7 GeV at W = 176 GeV and improves to 4.5 GeV at 
W = 225 GeV. 

The mass Mx of the dissociated photon system was reconstructed || by combining the 
information from the LUMI-e and the CAL: 



M x = VE* - P 2 « y/(E - P z ) ■ {E + P z ) = pE 1 ■ {E + P z ), (1) 

where E is the energy of the hadronic system observed in the CAL; the total momentum of the 
hadronic system, P, approximately equals the longitudinal component, Pz, as the transverse 



component generally is small in the case of photoproduction events. The following formula was 
used for the mass reconstruction: 



M x rec = 01 • /2£ 7 • (53 E t + £ Ei cos 0<) + o 2 . (2) 

y cond cond 

The quantities Ei and ^ denote the energy and the polar angle of CAL condensates, defined as 
groups of adjacent cells with total energy of at least 100 MeV, if all the cells belong to the EMC, 
or 200 MeV otherwise. These cuts reduce the effect of noise on the mass reconstruction. They 
were applied in addition to a noise suppression algorithm which discarded all EMC (HAC) 
cells with energy below 60 MeV (110 MeV); for isolated cells the thresholds were increased 
to 80 MeV (120 MeV). The coefficients a\ and a 2 correct for the effects of energy loss in the 
inactive material in front of the CAL and of energy deposits below the threshold. Their values, 
ai = 1.14 and a 2 = 1.2 GeV, were taken from ||. The masses in the range 4 < Mx < 40 GeV 
are reconstructed with a resolution aM x /M x ~ 0.8/wMx(GeV) and an offset smaller than 
0.5 GeV @. 

The variable t = (p — p') 2 , where p and p' are the incoming and the scattered proton four- 
momenta, respectively, can be evaluated as t « — (p|./xl)[1 + {M 2 /p\){l — x_l) 2 ]. Both pt and 
xl were measured with the LPS. For the data considered here, which have xl close to unity, the 
approximation t « —Pt/xl was used. Since, as mentioned earlier, the incoming proton beam 
has an intrinsic transverse momentum spread of a Px ~ 40 MeV and a PY ~ 90 MeV, which is 
much larger than the LPS resolution in transverse momentum, the measured value of t is given 
by the convolution of the true t distribution and the effect of the beam spread. Because of this 
we make a distinction between the true value of t and the measured value, t apparent = —p\jXL- 



3.3 Offline selection 

To select the final sample, the following conditions were imposed on the reconstructed data: 

• A scattered positron in the LUMI-e with energy in the range 12 < E' e < 18 GeV, corre- 
sponding to 176 < W < 225 GeV. 

• An interaction vertex reconstructed by the tracking detectors. 

• Mass of the dissociated photon system in the range 4 < Mxrec < 32 GeV. The lower 
limit eliminates the region dominated by resonance production; it also reflects the lower 
limit of Mx = 1.7 GeV for which Monte Carlo events were generated (cf. section f|). The 
upper limit is a consequence of the limit on xl (see below), since M x ~ W 2 (1 — xl). 



In addition, the detection of a high momentum proton in the LPS was required JTT| : 



One track in the LPS with xl > 0.97 was required. This is used to select diffractive 
events in which the proton remains intact. 



Protons with reconstructed trajectories closer than 0.5 mm to the wall of the beam-pipe, 
at any point between the vertex and the last station hit, were rejected. This eliminates any 
sensitivity of the acceptance to possible misalignments of the HERA beam-pipe elements. 
In addition badly reconstructed tracks are removed. 

The pt range was restricted to the interval 0.27 < px < 0.63 GeV, thereby removing 
regions where the acceptance of the LPS is very small or changes rapidly |[T] . 



After these selections, 641 events remained. 

3.4 Background 

The background contamination in the sample was mainly due to two sources. 

1. Some activity in the RCAL can accidentally overlap with the scattered positron of a 
bremsstrahlung event (ep — > e'jp) in the LUMI-e [|J. A large fraction of this back- 
ground can be identified since the bremsstrahlung photon is accepted by the LUMI-7. 
For bremsstrahlung, one has E' e + E y = E e , where E^ is the energy of the radiated pho- 
ton. For such events the energy deposits in the LUMI-e and the LUMI-7 calorimeters 
thus sum up to the positron beam energy. These events were removed. 

The unidentified events were statistically subtracted by including the identified back- 
ground events with negative weights in all the distributions, thereby compensating for 
the unidentified background events |26], ^7|. This subtraction was less than 



/o. 



A proton beam-halo track in the LPS can accidentally overlap with an event satisfying 
the trigger and the selection cuts applied to the variables measured with the central ZEUS 
detector (beam-halo event). The term beam-halo track refers to a proton with energy 
close to that of the beam originating from an interaction of a beam proton with the 
residual gas in the pipe or with the beam collimator jaws. Obviously, a beam-halo track 
is uncorrelated with the activity in the central ZEUS detector. For such a beam-halo 
event, energy and momentum conservation are not necessarily satisfied; in particular the 
quantity (E + Pz + 2P^ PS ) , where Pz PS is the Z component of the proton momentum 
measured in the LPS, may exceed the kinematic limit of 1640 GeV. The condition (E + 
Pz + 2P^ PS ) > 1655 GeV (thereby including the effects of resolution) identifies such 
events, which thus were rejected. 

In order to evaluate the residual contamination after all cuts, the distribution of 2P^ PS 
for identified beam-halo events was randomly mixed with the (E + Pz) distribution for 
all events, so as to create a distribution of (E + Pz + 2P^ PS ) for halo events. The 
observed (E + Pz + 2P^ PS ) distribution was then fitted as the sum of the diffractive 
Monte Carlo contribution (cf. section |j) and the beam-halo contribution just discussed. 
The relative normalisation of the two terms was left as a free parameter. For events with 
(E + Pz + 2Pz PS ) < 1655 GeV, the fraction of beam-halo events thus was estimated to 
be (6.3 ± 1.2)%. Here again the identified background events were included with negative 
weights in all the distributions in order to compensate for the unidentified beam-halo 
events. 



The number of events remaining after the background subtraction, i.e. after removing the 
identified bremsstrahlung and beam-halo events and after including the effect of the negative 
weights, was 515. The contribution from non-single diffractive dissociation processes, e.g. dou- 
ble diffractive dissociation, is expected to be of the order of a few per cent || and was not 
subtracted. 

Figure [I] shows the observed W, M x , xl and px distributions for the selected events af- 
ter background subtraction and including the correction for the effects of the trigger prescale 
factors. 



4 Monte Carlo simulation and acceptance determina- 
tion 



The reaction ep — > eXp was simulated using a Monte Carlo generator 28 based on a model 



calculation by Nikolaev and Zakharov |29| . The generated Mx distribution was reweighted with 
the sum of a pomeron-pomeron-pomeron [[Jjj contribution (da/dM x oc 1/M X ) and a pomeron- 
pomeron-reggeon contribution (da/dM x oc 1/M X ), so as to obtain a satisfactory agreement 
between data and Monte Carlo. As discussed in sect. |5|, however, the present results on the t 
distribution are largely independent of the details of the Mx spectrum simulation in the mass 
range considered in the analysis. 

All generated events were passed through the standard ZEUS detector simulation, based on 
the GEANT program [R(J, and through the trigger simulation package. The simulation also 
includes the geometry of the beam-pipe apertures, the HERA magnets and their fields. The 
spread of the interaction vertex position was also simulated and so were the proton beam angle 
with respect to the nominal direction and its dispersion at the interaction point. The simulated 
events were then passed through the same reconstruction and analysis programs as the data. 
In Fig. |I] the distributions for the reconstructed Monte Carlo events as a function of W, Mx, 
Xl and px are compared with those of the data. The distributions of the simulated events were 
normalised to the observed number of events, corrected for the effects of the prescale and of the 
background subtraction. The agreement between the data and the Monte Carlo distributions 
is satisfactory. 

The acceptance was computed as the ratio of the number of reconstructed Monte Carlo events 
in a bin of a given variable and the number of generated events in that bin. The acceptance thus 
includes the effects of the geometric acceptance of the apparatus, its efficiency and resolution, 
as well as the trigger and reconstruction efficiencies. 



5 Results and discussion 

The acceptance corrected t distribution, dN/d\t\, is shown in Fig. || It was obtained by 
correcting the measured t apparen t distribution bin by bin with the acceptance determined from 
the Monte Carlo simulation described above. 

The data were fitted with the function: 



%-*■*"> (3) 

where A is a constant. The resulting value of the t-slope is 

b = 6.8 ± 0.9 (stat.)il:? ( s Y st -) GeV" 2 , 

in the kinematic region 4 < M x < 32 GeV, 0.073 < \t\ < 0.40 GeV 2 and 176 < W < 225 GeV. 
In this region, the average value of W is 200 GeV and the average value of Mx is 11 GeV. The 
result of the fit is indicated by the continuous line on Fig. |^. 

The analysis was repeated for the two Mx ranges 4 < Mx < 8 GeV and 8 < Mx < 32 GeV. 
The results, b = 7.0 ± 1.3 GeV" 2 and b = 6.5 ± 1.3 GeV" 2 , indicate no variation with Mx 
within the present sensitivity. 

The quoted systematic uncertainty on b (Ab) was obtained by modifying the requirements 
and the analysis procedures as listed below: 

1. Sensitivity to the selection of the proton track (cf. fll||): 

• The sensitivity to the proton beam tilt with respect to the nominal was evaluated 
by systematically shifting px by 10 MeV. 

• The track selection requirements were tightened. 

• Events with p x PS > and with p x PS < were analysed separately, as a check of 
possible relative rotations of the LPS detector stations. 

• The data were divided into a "low acceptance" and a "large acceptance" sample 
depending on the position of the LPS stations. The latter varied slightly from run 
to run. 

The last three contributions dominate; by summing all four in quadrature, Ab = ±1.0 GeV" 2 
was obtained. 

2. Sensitivity to the other selection cuts and acceptance corrections: 

• The W range was restricted to 195 < W < 215 GeV, leading to b = 6.7± 1.4 GeV" 2 . 

• The Mx distribution in the Monte Carlo was varied between da/dM x oc (1/Mx) 1 ' 5 
and da/dM x oc (1/Mx) 3 . The corresponding variation of b was at most ±0.2 GeV" 2 . 

• The vertex requirement was removed or restricted to \Z vertex \ < 50 cm. The effect 
on b was at most ±0.4 GeV" 2 . 

Summing these contributions to Ab in quadrature yields Ab = ±0.5 GeV" 2 . 

3. Background corrections: the size of the beam-halo background was varied by two standard 
deviations, yielding negligible effects. The bremsstrahlung background correction was 
removed altogether, causing changes in the result smaller than 0.1 GeV" 2 . 

4. Evaluation of the t-slope: 



The t-slope was determined with an alternative method discussed in detail in refs. [11, 



31] • One can express the t apparen t distribution as a convolution of equation (|3|) and a 
two-dimensional Gaussian distribution representing the beam transverse momentum 
distribution. The slope parameter b can then be determined by fitting the convolu- 
tion of eq. (|) and the two-dimensional Gaussian to the acceptance corrected t appareni 
distribution. This method has the advantage that the data can be binned in t apparent 
for which the resolution is much better than for t, as discussed earlier. The value of 
b thus obtained was b = 7.3 ± 0.9 GeV -2 . 

• The t-slope was also obtained with a third method: the t apparent distributions in the 
data and in the Monte Carlo were compared and the Monte Carlo t distribution at 
the generator level was reweighted until the x 2 °f the comparison between data and 
Monte Carlo reached a minimum. The result thus found differed from the nominal 
one by less than 0.1 GeV -2 . 

• The sensitivity to the binning in t was studied by using an unbinned maximum 
likelihood method for the fit, which gave a result differing from the nominal one by 
less than 0.1 GeV 2 . 

The quadratic sum of these effects contributes Ab =i$|i GeV -2 . 

All contributions were summed in quadrature, yielding a total systematic error Ab =-\' 2 GeV -2 . 

Table 1 lists the present result together with those of the Fermilab photoproduction exper- 
iment E612 |] and that obtained by ZEUS for (Q 2 ) = 8 GeV 2 ||. The present result agrees 
within errors with both. This agreement suggests that at fixed W there is little dependence of 
the slope on Q 2 and that for real photons the W dependence is not strong. Note that the t 
range of our measurement is different from that of M, which is 0.02 < \t\ < 0.1 GeV 2 . A direct 



comparison should therefore be made with caution: for example, in elastic Tip scattering |32 
the t-slope measured in the range 0.1 < \t\ < 0.4 GeV 2 is about 1.2 GeV -2 lower than in the 
range covered by |l|; this difference is of the same size as the errors of our measurement. 

The weak Q 2 dependence of the t-slope in diffractive photon dissociation may be contrasted 
with the change with Q 2 observed for elastic p° meson production [9-13], where b decreases by 
approximately 3-5 GeV -2 when going from Q 2 rj to Q 2 ~ 10 GeV 2 . If factorisation of the 
diffractive vertices |32| is assumed, the amplitudes for the reactions jp — > Xp and 7p — > p°p 
are proportional to the products of vertex functions G^x(Q 2 , t) ■ G pp {t) and G lp o(Q 2 , t) ■ G pp (t), 
respectively. In this framework, the t-slope includes the sum of the contributions from the j-X 
or 7-p vertex and from the p-p vertex. The comparison of the 'jp — > Xp and 7p — > p°p slopes 
indicates that the vertex function G.yx{Q 2 ,t) has a weaker dependence on Q 2 than G lp o(Q 2 ,t). 
A rapid Q 2 dependence of G lp o(Q 2 ,t) is expected in pQCD inspired models of elastic vector 
meson production |33[], reflecting the decrease with Q 2 of the transverse size of the quark- 



antiquark pair into which the photon fluctuates before interacting with the proton. A weak Q 2 
dependence of the function G^xiQ 2 ^) is also expected in the framework of various models of 
diffractive dissociation of photons 
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<Q 2 )/GcV 2 
W range/GeV 
Mx range/GeV 
\t\ range/GeV 2 
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8 

50-270 
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6/GeV- 2 


6.8 ±0.9 (stat.) t\ 2 x (syst.) 


4.2 ±1.4, 6.3 ±1.3, 5.1 ±1.3 


7.2 ±1.1 (stat.) to.l (syst.) 
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is listed together with those from ref. QXJ] for real photons and that of ref. || for virtual photons. 
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Figure 1: Distributions as a function of (a) W, (b) Mx, (c) £l and (d) pr for the data (points) 
and the Monte Carlo events (histogram). The error bars indicate the statistical uncertainty 
of the data. The number of data events is corrected for the trigger prescale factors and the 
background. 
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Figure 2: Differential distribution dN/d\t\ for photon diffractive dissociation, '-/p — > Xp, in the 
kinematic region 176 < W < 225 GeV and 4 < Mx < 32 GeV. The vertical bars indicate the 
size of the statistical uncertainties. The line is the result of the fit described in the text. The 
scale on the vertical axis is arbitrary. 
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